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Jet Aero acoustics: Noise Generation Mechanism and Prediction 

Progress Report 

This report covers the third year research effort of the project. The research work 
focussed on the fine scale mixing noise of both subsonic and supersonic jets and the 
effects of nozzle geometry and tabs on subsonic jet noise. Technical results of the 
research effort are reported in the following publications. 

1. Tam, C.K.W. and Auriault, L. "Jet mixing noise from fine scale turbulence", 
AIAA Paper 98-2354, to appear in the AIAA J. Feb. 1999. 

2. Tam, C.K.W. and Zaman, K.B.M., "Subsonic jet noise from non-axisymmetric 
and tabbed nozzles", AIAA Paper 99-0077, 1999. 

(Copies of the above are attached to the end of this report. 

In publication 1, a new semi-empirical theory of jet mixing noise from fine scale 
turbulence is developed. By an analogy to gas kinetic theory, it is shown that the 
source of noise is related to the time fluctuations of the turbulence kinetic theory. 
On starting with the Reynolds Averaged Navier-Stokes equations, a formula for the 
radiated noise is derived. An empirical model of the space-time correlation 
function of the turbulence kinetic energy is adopted. The form of the model is in 
good agreement with the space-time two-point velocity correlation function 
measured by Davies and coworkers. The parameters of the correlation are related to 
the parameters of the k — e turbulence model. Thus the theory is self-contained. 
Extensive comparisons between the computed noise spectrum of the theory and 
experimental measured have been carried out. The parameters include jet Mach 
number from 0.3 to 2.0 and temperature ratio from 1.0 to 4.8. Excellent agreements 
are found in the spectrum shape, noise intensity and directivity. It is envisaged that 
the theory would supercede all semi-empirical and totally empirical jet noise 
prediction methods in current use. 



Publication 2 uses the two similarity spectra discovered by Tam, Golebiowski and 
Seiner to demonstrate that the jet noise fields from circular, elliptic or rectangular 
nozzles are nearly identical. Nozzle shape modification is not an effective way to 
suppress jet noise. The effects of inserting tabs into a subsonic jet is also 
investigated. A jetlets model is developed to correlate the noise from a jet with N 
tabs to that of a single jet. The jetlets model predicts that the noise spectrum is 
shifted to a higher frequency by a factor of N in . This result is in excellent agreement 
with experiments. This theory is the only theory on the noise of jets from tabbed 
nozzles in the literature. 
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Abstract 

It is known that turbulent mixing noise from 
high-speed jets consists of two components. They 
are the noise from large turbulent structures in the 
form of Mach wave radiation and the less directional 
fine scale turbulence noise. The Mach wave radia- 
tion dominates in the downstream direction. The 
fine scale turbulence noise dominates in the side- 
line and upstream directions. In this paper, a semi- 
empirical theory is developed for the prediction of 
the spectrum, intensity and directivity of the fine 
scale turbulence noise. The prediction method is 
self-contained. The turbulence information is sup- 
plied by the k — e turbulence model. The theory 
contains three empirical constants beyond those of 
the k - e model. These constants are determined 
by best fit of the calculated noise spectra to exper- 
imental measurements. Extensive comparisons be- 
tween calculated and measured noise spectra over 
a wide range of directions of radiation, jet veloci- 
ties and temperatures have been carried out. Excel- 
lent agreements are found. It is believed that the 
present theory offers significant improvements over 
current empirical or semi-empirical jet noise predic- 
tion methods in use. There is no first principle jet 
noise theory at the present time. 

1. Introduction 

Since the pioneering work of Lighthill 1,2 , there 
have been numerous attempts to develop a jet noise 
prediction theory. In the literature, one finds many 
proposed theories and semi-empirical theories. The 
main difficulty in predicting jet noise is our lack of 
understanding of the turbulence in jet flows. This is 
true even today. However, many of the fundamen- 
tal physics of jet flows that affect the propagation 
and radiation of jet noise were recognized very early. 
They have been incorporated into some of the the- 
ories. These effects include mean flow refraction 3,4 
and source motion 5 . An excellent review on jet noise 
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theories up to the late 1980 ’s can be found in Ref. 
[6]. Ref. [7] to [9] provide reviews on the same sub- 
ject. However, they include some of the more recent 
works and the emphasis and perspectives are quite 
different. 

Recently, Tam, Golebiowski and Seiner 10 sug- 
gested that, because there is no intrinsic charac- 
teristic length and time scales in the mixing layer 
of a high Reynolds number jet (up to the end of 
the core region), not only the mean flow and the 
turbulence statistics must exhibit self-similarity, the 
same must be true for the radiated noise. By ex- 
amining the entire data bank (1900 spectra in all) 
of the Jet Noise Laboratory of the NASA Langley 
Research Center, they found that turbulent mix- 
ing noise of high-speed jets consisted of two distinct 
components (see also Ref. [9]). Each component ex- 
hibits self-similarity of its own. One component ra- 
diates principally in the downstream direction. This 
is consistent with Mach wave radiation from the 
large turbulence structures/instability waves of the 
jet flow 11 " 13 . The other component that has a rela- 
tively uniform directivity is dominant in the sideline 
and upstream directions. These characteristics sug- 
gest that it is the noise from the fine scale turbulence 
of the jet flow. Tam et al. 10 succeeded in identifying 
two similarity spectra from the data. They demon- 
strated that one of the spectra they found fitted the 
noise from the large turbulence structures and the 
other the noise from the fine scale turbulence regard- 
less of the jet velocity, temperature and direction of 
radiation. More recently, Tam 14 showed that even 
the noise spectra of non-axisymmetric jets including 
jets from rectangular, elliptic, plug and suppressor 
nozzles fitted the same two empirically found simi- 
larity spectra. This indicates that the noise sources 
of these jets are similar to those of the circular jet. 
Their noise is also made up of two components. 

The primary objective of this work is to develop 
a semi-empirical theory for the prediction of the fine 
scale turbulence noise from high-speed jets. The 
present theory is self-contained. The turbulence in- 
formation of the jet flow needed for noise prediction 
are supplied by the k— e turbulence model. Since the 
k — e model is a semi-empirical model, the present 
theory is also semi-empirical. As will become clear 
later, the noise prediction formula developed here 
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contains three empirical constants. These constants 
are determined by best fit to the noise data. Once 
the constants are decided, the formula can provide 
accurate noise spectrum and directivity predictions 
over the jet velocity ratio range of — 0.4 to 3.0 
( uj and aoo are the fully expanded jet velocity and 
ambient sound speed) and temperature ratio range 
of If*- = 1.0 to 5.0 (T r and Too are the jet reservoir 
and ambient temperature). This parametric range 
brackets all known jet noise experiments and com- 
mercial jet engine operating conditions. 

The use of the k — e turbulence model to predict 
jet noise is not new. Khavaran, Krejsa and Kim 15,16 , 
Bailly and coworkers 17 ” 20 employed the k—e model 
to provide turbulence statistics to their chosen ver- 
sion of the acoustic analogy theory for supersonic jet 
noise prediction. Khavaran et ai tuned their empir- 
ical constants to fit the measured jet noise data at 
Mach 1.4. However, even at this Mach number, the 
calculated spectra do not fit the measured data well 
at 90 degrees and in the forward direction for which 
the fine scale turbulence noise dominates. Bailly and 
coworkers were more interested in using their model 
to calculate Mach wave radiation from supersonic 
jets. There is fair agreement between their com- 
puted and measured directivity data but the calcu- 
lated spectral distributions do not match well with 
experiments. 

To develop a fine scale turbulence noise theory 
requires many steps. We will first turn to the ques- 
tion of how fine scale turbulence in a jet generates 
sound. For this purpose, we will use the gas kinetic 
theory analogy. The use of gas kinetic theory to 
elucidate ideas in turbulence is well established 21 . 
Consider the motion of gas molecules in a moving 
frame fixed to the mean velocity v as shown in fig- 
ure 1. Suppose m is the mass of a molecule, n is 
the number density. Due to the random motion of 
the gas molecules, the gas exerts a pressure, p, on 
its surroundings. It is a simple matter to show, fol- 
lowing standard kinetic theory of gases 22,23 , that p 
is given by, 

p= ^mn(v-v) = ^p(v 2 ) (1) 

where v is the random molecular velocity, p is the 
density of the gas and ( ) is the ensemble average. 

Now we will regard fine scale turbulence as small 
blobs of fluid moving randomly as shown in figure 2. 
Again, let v be the random velocity of the fine scale 
turbulence measured in the mean flow moving frame. 
By analogy to the gas molecules of figure 1, the fine 
scale turbulence effectively exerts a pressure, p tur b» 


on its surroundings. This pressure, following (1), is 
equal to 

1 2 

Pturb = 9* = 3 p(*> 2 ) = ^pk, (2) 

where k 3 = \{v 2 ) is the kinetic energy of the fine 
scale turbulence per unit mass. The pressure given 
by (2) is a macroscopic quantity valid for length scale 
larger than the size of an individual blob of turbu- 
lent fluid. This pressure must be balanced by pres- 
sure and momentum flux of the surrounding fluid. 
If this pressure fluctuates in time, it will inevitably 
give rise to compressions and rarefactions in the fluid 
medium. This results in acoustic disturbances. Fol- 
lowing this reasoning, one would expect the source 
of fine scale turbulence noise to be equal to the time 
rate of change of Pturb or Qs in the moving frame 
of the fluid. In other words, the convective 

derivative of q s , is the noise source term in a turbu- 
lent flow. 

In section 2, the equations governing the gener- 
ation and propagation of acoustic disturbances in a 
jet flow (generated by fine scale turbulence) are for- 
mally derived. Here we will begin with the Reynolds 
Averaged Navier-Stokes equations (RANS). The re- 
lationship between the radiated sound field and the 
source of fine scale turbulence noise will be estab- 
lished. In section 3, a semi-empirical model of the 
noise source space-time correlation function is de- 
veloped. The parameters of the correlation function 
are then related to the length and time scales of the 
k — e turbulence model. By means of this source cor- 
relation function, a formula for the spectrum of the 
radiated noise is derived. Extensive comparisons be- 
tween the calculated results of this formula and ex- 
perimental measurements have been carried out. A 
sample of these comparisons are reported in section 
4. Good agreements are found over both subsonic 
and supersonic Mach number. Temperature depen- 
dence over a wide range of jet to ambient tempera- 
ture ratios is also correctly predicted. 

2. Formulation 

Consider the flow of a turbulent compressible 
jet. A convenient way to account for the effect of 
compressibility is to use Favre average variables 24 . 
For the present purpose, we may regard the average 
as a volume average 25 with a sharp Fourier cutoff 
filter; the volume being smaller than the large tur- 
bulence structures of the jet flow but larger than 
the fine scale turbulene. We will decompose a flow 
variable into two parts, 

/ = /+/" ( 3 ) 


2 



where / and /" are the Favre average and Favre 
fluctuating components. On starting from the equa- 
tions of motion, it is straightforward to derive the 
familiar Reynolds averaged Navier-Stokes equations 
(RANS); e.g., see Ref. [26]. 


P 


fdui du { 
\ dt Uj dxj 


) 


dp d 
dx{ dxj 


(-pu'/tiV + Tfjj 


(4) 


where an overbar denotes a Reyn olds averaged quan- 
tity. The terms T{j and —pu'/uj on the right side of 
(4) are the viscous and Reynolds stresses, respec- 
tively. On following the generally accepted RANS 
approach, we will adopt the Boussinesq eddy viscos- 
ity model for the Reynolds stresses; i.e., 


-pu'/u'J = 2 Pt (sij - - | pk s Sij 

where pt is the eddy viscosity and 

l fdui duA A _du k 
bij ~ 2\dxj + dxj)’ 


dx k 


-pu'/u'S = 2 pk } . 


(5) 

( 6 ) 
(7) 


Of importance is that k, as defined by (7) is the 
kinetic energy of the fine scale turbulence per unit 
mass. 

As discussed in the previous section, k s in (5) 
effectively represents a pressure field exerted by the 
fine scale turbulence on the surrounding fluid. The 
time fluctuation of this pressure field causes com- 
pressions and rarefactions in the fluid medium. In 
this way, sound is generated by the small scale tur- 
bulence. Inside the jet flow, the acoustic field is a 
very small part of the unsteady fluctuations. It is 
sufficient to consider the acoustic field generated by 
the time dependent part of k Si denoted by k 9 , to 
be given by the linearized form of equations (4) and 
(5). The linearization is to be performed over the 
mean flow of the jet. In addition, both molecular 
and eddy viscosity terms will be ignored as they will 
have only a relatively small effect on the acoustic 
disturbances. The linearized form of (4) with only 
k s terms retained on the right side, after applying 
the locally parallel mean flow approximation, is 


_ [ dui dui duj 

P\ — + uj— + «i — “ 


dt 


dxj 


dxj 
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dp 

dxi 


dg, 

dxi 


where 


2 -t 

<h = 3 pk a 


(8) 


(9) 


is the unsteady pressure exerted by the fine scale 
turbulence. In (8), tq and p are the mean veloc- 
ity and density of the jet and iq, p are the acoustic 
field variables associated with k s . We will now sup- 
plement equation (8) with the linearized continuity 
and energy equation (ignoring viscous and thermal 
dissipations) to form a closed system of equations 
for the acoustic field. With respect to a cylindri- 
cal coordinate system (r, <f>,x) centered at the nozzle 
exit (with the x-axis coinciding with the jet axis), 
the full set of governing equations are, 


du __ du du 
-^7 + ti ■=- + v— 


dt 


dx 


dv _dv 1 
dt +U !h 
dw _ dw 
~dt +U !h 


dp _ dp _ 


dr 

dr 
1 dp 

+ "77 = 

r d<p 
*1 d(vr) 
r dr 


_ dg s 
dx dx 
dg, 
dr 

- 1 dg 3 

r d<j> 

1 dw du 
^ r d<f> + dx 


(10a) 

(106) 

(10c) 


= 0. 


(10cO 


In these equations (u,v,u>) are the velocity compo- 
nents in the (x,r,<£) directions. Also, in accordance 
of the locally parallel flow approximation, p, u are 
regarded as functions of r alone and p = Pco is con- 
stant. 


2.1. The Acoustic Field 

Equations (10a) to (lOd), except for the noise 
source terms on the right, are identical to the Lil- 
ley’s equation 4,6 . As is well known, the Lilley’s ap- 
proach was designed to account for the mean flow 
refraction effect. This effect is especially significant 
for the high frequency part of the noise spectrum. 
This is confirmed in the numerical results reported 
in section 4. 

A formal solution of (10) directly relating the 
radiated acoustic pressure field to the fine scale tur- 
bulence intensity q s can be found by using the space- 
time Green’s function of these equations. Since q 3 
appears on the right side of equations (10a), (10b) 
and (10c), three Green’s functions are necessary. Let 
(um^n^njPn) with n = 1, 2 and 3 be the Green’s 
functions corresponding to a source on the right 
side of equations (10a), (10b) or (10c), respectively. 
These Green’s functions satisfy the following nonho- 
mogeneous equations. 


du n _ du n du 

-m +u + 


+ 


dpn 

dx 


= S(x-xi)S(t - U)6 nl 


( 11 ) 
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( 12 ) 


_ \dv n _ 3u n ] dp n 

/> ["aT + u 

= 8{x-xi)6(t —ti)5„2 

_ [ dw„ _ dw n ] 1 dp n 

P [~dr + U ~d^\ + r ~d* 


= 8(x-xi)6(t-ti)8 n3 
1 d{v n r) 
r dr 

= 0 


dp n . - dpn . - 

-m + “ a7 + ip 


(13) 


1 dw n du n 

r d<j> dx 


(14) 


where n — 1,2,3. Here J( ) is the Dirac delta func- 
tion and 6 nm is the Kronecker delta. The Green’s 
functions depend on four variables. They are the 
field or observer coordinates, x, the source coordi- 
nates, Xi, the field or observation time, t , and the 
source time, t\. To avoid confusion, whenever it is 
necessary to display the arguments of the Green’s 
function, they will always be in exactly this order; 
e.g., p n (x,Xi,Mi). 

The space-time Green’s function is related to the 
time harmonic Green’s function (u n ,v n ,w n p n ) by 
the Fourier inverse transform. For example, we have 


functions are related to the adjoint Green’s function 

by 27 , 


Pl(x, Xi, w) = ii a (x 1 ,x,w) 

P 2 (x,x i,w) = v a (xi,x,u ;) (20) 

P3( x > x i> w ) = WafXi.X.w). 


It is worthwhile to point out that the arguments on 
the two sides of equation (20) satisfy the reciprocity 
relation. That is, the source point and the field point 
of the adjoint are interchanged. It is clear from (20) 
that there is another important advantage in using 
the adjoint Green’s function. If one is interested in 
the far acoustic field, namely, the pressure field, it is 
only necessary to solve the adjoint equations (16) to 
(19) once to obtain (u a , ^a> ^a)- On the other hand, 
to find fa, p 2 and p 3 directly one has to solve the 
original time harmonic Green’s function equations 
three times. 

By means of the adjoint Green’s function, the 
pressure field generated by the source terms on the 
right side of (10a) to (10d) is formally given by, 


oo 

Pn(x,Xi,Mi) = J p n (x,Xi,w)e“’ u ' (t " tl) du; (15) 

— OO 


The equations governing the time harmonic Green’s 
functions are simply the Fourier transforms of (11) 
to (14). Recently the present authors 27 have shown 
that there is great computational advantage to use 
the adjoint time harmonic Green’s function instead 
of the time harmonic Green’s functions. On fol- 
lowing the work of Ref. [27], it is easy to find 
that the adjoint time harmonic Green’s function 
(u a , v a , Way Pa) (subscript a indicates the adjoint) 
satisfies the following equations. 


_ du a 
dx 


~ p j tUJUa 4- u 
_ [. _ dv a du 

~ f l"”* + *“‘ 


-dp a r 


(16) 


- 0 < 17 > 


- p 


. -SlOal 

lUJW a + u 


dx 


_ TP_dpa_ _ 0 
r d<f> 


(18) 


_ dp a 


"1 d[v a r) 1 dw a dug 
r dr + r d<j> dx _ 


= s i(x-x 1 ) 


(19) 


It is straightforward to show that the pressure 
fields of the three original time harmonic Green’s 


00 

P(X,<) = -////// [..(x„x, U ) 

— CO 

.^(xi.t,) 


dqAxuti) 

dx i 


+ Vo(Xl,X,w)- 


+ 


dr i 

^a(xi,x,u;) dqsjxiJi) 


r d<f> 

e -i u (t- tl ) d ijdtidxi. 

( 21 ) 

Let Vi = (afr e x + 9|7 e y + a7: ez ) be the 6 radient 
with respect to the xi coordinates. (21) may be 
rewritten in the form, 


p(M) = - 


//////' [v '' (v * (Xl 
— OO 




-qs(*i,ti) Vi •v a (xi,x,w)] (22) 

. e -Mt-ti)d u dhdxi 


where v a = (u a> v a , w a ). This formula can be simpli- 
fied in two ways. First, we can apply the divergence 
theorem to the first term of the integrand. This con- 
verts the volume integral dxi to a surface integral 
outside the jet flow. But q s is zero outside the jet. 
Thus, there is no contribution from the first term. 
Now for v a (xi,x, w), the source is outside the jet flow 
if x is in the far field. Thus by making use of (19), 
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(22) becomes, 


p(m) =_ ////// 9 *( x i’*i) (* w+u £r) 

x outside the jet J J J J L \ u 1 / 

dhdxx 

-llfllb >l(s-s:) 


dtxdxi . 

Finally, by simple integration by parts, we find the 
desirable formula for p(x,<), 


00 00 

p(x,t) =11111 / p a (x I,x,w) 

itside the iet J J J J J J 


x outside the jet 


~oo L-00 


(23) 


. e - Mt-ttifo D( l^\ lhldt l dx l 
Dt 1 


where + u is the convective derivative 

following the mean flow. 

Equation (23) is the main result of this section. 
Obviously, it is open to the interpretation that noise 
from fine scale turbulence is generated by the time 
rate of change of the turbulence kinetic energy or 
pressure in the moving frame following the convec- 
tion of the fine scale turbulence by the mean flow. 
The term inside the square bracket is the adjoint 
Green’s function with the source located at x and 
the field point at Xi according to the reciprocity re- 
lation. 


2.2. The Spectral Density of the Radiated 
Sound 


By means of (23), the autocorrelation function 
for a point in the acoustic far field may be formed; 
i.e., 


00 

<p(x,f)p(x,f + r)) = JJ-JJp a (Xi,X,U>i) 


, „ / Dq^Xutx) Dq,{x 2 ,t2)\ 

— ; (24) 


In (24), ( ) is the ensemble average. 

The spectral density of the radiated sound, 
S(x,w) is the Fourier transform of the autocorre- 
lation function or 


00 

S(x,u)= i- J {p(x,t)p(x,t + r))e iuT dr. (25) 

— OO 


OQ 

By using f exp[i(w — = 2ttS(uj — W2), it 

—00 

is straightforward to derive from (24), the following 
formula for S(x,w). 


00 

S(x,w) = JJ-JJp a (xi,X,Ui)p a (X 2 ,X,U 2 ) 

— OO 

/ Dq,(*l,tl) £ ) ? J (X2,<2) \ 

\ Dt\ Di 2 J 

• 5(u — W2)duJiduj2dtidt2dxi dx 0 . 


(26) 


3. Model Space-Time Correlation Function 

To proceed further, we need a mathematical rep- 
resentation of the noise source space-time correla- 
tion function This correla- 

tion function has never been measured before. How- 
ever, two-point space-time correlation of the fluctu- 
ating axial velocity component in jets has been mea- 
sured and studied by Davies, Fisher and Barrett 28 
and Chu 29 . We believe that the mathematical form 
of {7^ 7^) should be similar to that of the mea- 
sured two-point space-time correlation function of 
the axial velocity component. Here we propose to 
adopt a model space-time correlation function char- 
acterized by three parameters. We will show that 
with an appropriate choice of the parameters, the 
model function fits the measured function of Davies 
et ai (see figure 3) well. 

Let £ = xi - x 2 , T) = 2/1 - y 2 , C = *1 - *2, 
r — t x — t 2 . We will consider the following model 
function 

/ £M X l.*l) ^( x 2,<2 ) \ _ ^ 

\ Dti Dt 7 } c 2 r 2 (27) 

In (27) u is the mean flow velocity or the transport 
velocity of the fine scale turbulence. The three pa- 
rameters of the model are t 9 , r 9 and q 2 S . t 9 is the 
characteristic size of the fine scale turbulence in the 
moving frame of the mean flow. r 9 is the character- 
istic decay time. ur 9 is the decay distance. In the 
limit £ , *7, C T 0> (27) becomes 



Thus q s is a measure of the RMS value of the fluctu- 
ating kinetic energy of the fine scale turbulence and 
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ct 3 represents a typical time scale of the fluctuation. 
The coefficient c is expected to be less than 1.0. This 
indicates that the fluctuating time is shorter than 
the turbulence decay time. 

Figure 3 shows a plot of the correlation func- 
tion (27) with l s = 0.1758 inches, r s = 447.4/1 sec. 
and u = 3.515 x 10“ 3 inch //z sec. As can be seen, 
the model function is in reasonably good agreement 
with the measurements of Davies et al. 28 . This re- 
sult assures us that the model space-time correlation 
function does have the right functional character- 
istics. We will use this model correlation function 
throughout this work. 

Model correlation function (27) has three pa- 
rameters, namely, £ s , r s and q 3 . They are the char- 
acteristic parameters of the fine scale turbulence of 
the jet flow. Here we propose to obtain the values 
of these parameters through the k — e turbulence 
model. 

It is generally known that the standard k — e 
turbulence model does poorly in predicting the mean 
flow of jets 30 . Thies and Tam 31 recognized the prob- 
lem of applying the standard k — e model to jet flows. 
They proposed to modify the k—e model coefficients 
and showed convincingly by comparison with a large 
set of jet flow data over the jet Mach number range 
of 0.4 to 2.0 and temperature ratio of 1.0 (cold jet) 
to 4.0 that their modified k — e model is reliable and 
accurate. In this work, the modified k — e model of 
Ref. [31] is used. 

The k — e model provides only two pieces of in- 
formation about the turbulence of the jet flow. They 
are k , the averaged turbulence kinetic energy, and e, 
the dissipation rate. But with k and e known, it is 
possible to form a length i , characterizing the size 
of the small scale turbulence and a decay time, r, of 
the turbulence as follows, 


t = 




(29) 


The parameters £ and r of (29) are directly rele- 
vant to the parameters £ s and r, of model two-point 
space-time correlation function (27). To establish 
their relationship, we performed a k — e model cal- 
culation for the jet flow of the Davies et al. ex- 
periment. At the location of the measured correla- 
tion function shown in figure 3, the value of t and r 
obtained were very close to 0.176 inches and 447 /z 
sec., respectively. These are the values of t s and r s 
used to match correlation function (27) to the mea- 
surements. We, therefore, believe that £ s and r, of 
our model correlation function is directly related to 
t and r of the k — e model. However, the k — e 


model includes contributions from the large turbu- 
lence structures whereas £ s and r s are those of the 
fine scale turbulence alone. On accounting for this 
difference, we propose to let, 

„ , jt§ k 

£ s — C t i — c t — , r s — c T r — c T - (30) 

where c* and c T (c/,c T < 1.0) are constants to be 
determined empirically. 

In the above, we have noted that q 3 is a mea- 
sure of the intensity of the fluctuation of the kinetic 
energy of the fine scale turbulence. We, therefore, 
expect it to be proportional to q. (? = §pfc.) We 
will let 

4 = (31) 

c £ 

where A is the third empirical constant. 

With (30) and (31), the two-point space-time 
correlation function (27) is known save for the three 
constants c*, c T and A. These constants are to be 
determined by best fit of the calculated noise spectra 
to experimental measurements. 


4. A Formula for the Noise Spectrum 

On substitution of (27) into (26), the noise spec- 
trum in the far field is given by 


5(x,w) 



.X,Wl)pa(X2,X,UI 2 ) 


-2 

Qs 

c 2 rp 


- 1 -L ^ — t 3 )) 2 + (yi-y 2 ) 3 + (2i -*a)*l 

• e ’ • 


g — 1(0/1 +a> 3 )fc*f iwi t 


(32) 


• S[uj — U2)dtidt2du>2duJ2d Xi d X 2 


We will now show that most of the above integrals 
can be evaluated analytically. The entire expression 
can be reduced to a single volume integral over the 
jet flow. 

The first step is to integrate dt\ . The integration 
can best be carried out by first making a change of 
variable to s where 

Next, the <2 integration can be performed resulting 
in 27 tJ(cji H- u/ 2 ). With S(u 1 + u/ 2 ) and S(uj - 0 ^ 2 ) in 
the integrand, the u 1 and u? 2 integrals can be readily 
evaluated giving 


5(x,w) = 2 tt 


(n)*M£ 


u 


•Pa(xi,x,— u/)p a (x 2 ,x,w)e 

• t 


yif [x L -* 2 l 
u 3 (< In 2) “ r g 


»w(x! — X 2 ) 

“ dxid%2- 


(33) 
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In (33), the integrand is nearly zero unless x x is 
close to X 2 - Since x is in the far field, the adjoint 
Green’s function p a (xi,x, -u;) and p a (x 2 ,x, -w) dif- 
fer primarily by a phase factor. The phase difference 
can be assessed by ray acoustics. Figure 4 shows the 
rays from a far away point x to Xi and X 2 . The rays 
are essentially parallel except after entering the jet 
flow. We will now use a spherical polar coordinate 
system centered at the nozzle exit with the polar axis 
coinciding with the z-axis, or the jet axis. Let the 
spherical coordinates of x be (/?, 0, <f>) where 0 is the 
polar angle. The paths for the rays to Xi and X 2 differ 
by a length nearly equal to AC = (xi — X 2 ) cos 0 (see 
figure 4). The phase difference is equal to wavenum- 
ber times the difference in path lengths. Thus, the 
two adjoint Green’s functions differ by a phase factor 
of exp[i“- cos 0(xi — X 2 )] where aoo is the ambient 
sound speed. 

We will now adopt the approximation 

Pa(xi,x,-w) ~Po(X 2 ,X,-w)e'^ COS ® (l,_l3) . (34) 

(34) may also be obtained by invoking a locally 
parallel flow approximation. On substitution of 
(34) into (33) and on noting that Pa(*2, x > “ a? ) = 
p*(x 2 ,x,u;) where * indicates the complex conjugate, 
the integrals over Xi can be easily evaluated. This 
gives, 


S(x, w) = 47 t 



|Pa(x 2 ,x,cj)| 2 e_ 

[1 + w 2 r 2 (l — j 4 - cos 0) 2 ] 2 


(35) 


The spectral density S(x,u;) has the dimen- 
sions of (pressure) 2 time. On converting to dB per 
Strouhal number (^ i ; uj and Dj are the fully ex- 
panded jet velocity and diameter) the spectral den- 
sity of the sound pressure field at (/?,©,<£) is, 


S = 10 l °S 


47r5(x, w) 
2 Dj— 

^ref. uj 


(36) 


where p re f is the reference pressure of the decibel 
scale. 


5. Comparisons with Experiments 

Before we compare the numerical results of for- 
mulas (36) and (35) with experimental measure- 
ments, it is believed worthwhile to describe briefly 


how the various quantities and integrals of these for- 
mulas are determined and evaluated. Essentially the 
computations are carried out in three steps. First, 
with a prescribed jet operating condition, the mean 
flow and the values of k and e are calculated follow- 
ing Ref. [31]. Once k and e are found, £ S1 r s and ^ 
are determined by (30) and (31). The second step is 
to make use of the calculated mean flow to find the 
adjoint Green’s function as in Ref. [27]. Finally, the 
volume integral is evaluated by dividing the jet axi- 
ally into slices of 0.5 diameter. Azimuthally, the ad- 
joint Green’s function is expressed in Fourier modes. 
They are summed over. In the radial direction, a 
very fine mesh is used in the volume integration. 

We will compare formula (36) with four sets of 
data. They are the data of the Jet Noise Laboratory 
of the NASA Langley Research Center measured by 
Seiner 10 and his coworkers, the Norum & Brown 32 , 
Tanna et al 33 and Ahuja 34 data. The NASA Lan- 
gley data are from supersonic jets at temperature 
ratio of 1.0 to 4.9. The Norum & Brown data in- 
clude both subsonic and supersonic cold jets. The 
Tanna et al. data are from subsonic and supersonic 
jets. Only hot jet data at supersonic Mach num- 
ber are considered of good quality to be included 
for comparison. The Ahuja data are all at subsonic 
Mach number at room temperature. There is overall 
consistency among the four sets of data. However, 
for hot jets at supersonic speed, the NASA data is 
approximately 1 to 2 dB higher than the Tanna et 
al. measurements. Generally speaking, we believe 
that there is an uncertainty in the data of the order 
of 1 to 2 dB. 

Before formulas (35) and (36) can be used, the 
three empirical constants c*, c T and A have to be 
determined. After some trial and error to obtain 
best fit to the data, the following values are selected, 

ct = 0.256, c T = 0.233 and A — 0.755. (37) 

These are the values used throughout the rest of the 
paper. 

Figures 5a, b, c show comparisons of the calcu- 
lated noise spectra and the experimental measure- 
ments of Ref. [10] and [32] at = 1.0. The 
spectra are scaled to a distance of 100 Dj. Shown 
in these figures are data from jets at fully expanded 
jet Mach number Mj = 2.0, 1.49, 0.9, 0.7, 0.5 and 
0.3. Figure 5a is for noise radiation at 0 = 60 
deg. As can be seen, there are excellent agreements 
between the calculated spectra and measurements 
except at very high frequency. At very high fre- 
quency, the measurements are less reliable as they 
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are close to the high frequency cut-off of the mi- 
crophones. Figure 5b shows similar comparisons at 
0 = 90 deg. The agreements are again excellent ex- 
cept for the Mj — 1.49 case. The two broadband 
peaks on the high frequency side of the spectrum 
are broadband shock associated noise 35,7,8 . They 
should be ignored in the comparison. In Ref. [10], 
this was noticed when comparing the spectrum of 
a perfectly expanded jet against that of an imper- 
fectly expanded jet at identical operating conditions. 
Figure 5c shows the comparisons at 0 = 120 deg. 
Again, there are good agreements. In the forward di- 
rection, most supersonic jet spectra, invariably, con- 
tain broadband shock noise. The difference between 
the calculated and measured spectra at Mj — 2.0 
and 1.49 is mainly due to the presence of broadband 
shock associated noise. 

Figure 6 shows the comparisons of hot jet spec- 
tra. The jet Mach number and temperature ratio 
are fixed at 2.0 and 1.8, respectively. The spectra* 
in three directions are provided. There are good 
agreements between the calculated spectra and the 
experimental measurements. This is true both in the 
peak level and spectral shape. The slight difference 
for the case at 0 = 117 deg. is again due to the 
presence of broadband shock associated noise in the 
forward direction. 

To test whether formula (36) can predict accu- 
rately the change in noise spectrum with jet tem- 
perature, the noise data of Ref. [10] at Mach 2.0 in 
the © = 90 deg. direction are used for comparison. 
Figure 7 shows the calculated and measured spectra 
at four jet temperature ratios ranging from 1.0 to 
2.72. As can readily be seen, there is again good 
agreement in each case. 

To provide a global comparison of the calculated 
results of formula (36) and experimental measure- 
ments over a broad range of jet velocity and temper- 
ature, the peak value of S(R } 0, <t>\ J ^ L ) at 0 = 90 
deg. are plotted in figures 8a and 8b as a function 
of dLL, it was found in Ref. [10], the data in such a 
plot tended to align themselves along nearly parallel 
lines according to the jet temperature ratio Fig- 
ure 8a shows the comparison between the calculated 
results of formula (36) and the data from Ref. [10] 
and [32]. On considering that there is a fair amount 
of scattering in the experimental data, the overall 
agreement is excellent. Here the jet temperature ra- 
tio varies over a range of 5. This brackets all present 
day commercial jet engine operating conditions. The 
calculated noise level and measurements are in close 
agreement at nearly every temperature ratio. Fig- 
ure 8b shows the comparison between the calculated 


results and the measurements of Tanna et a/. 33 and 
Ahuja 34 . There is a fair amount of scattering of the 
subsonic jet data. But they cluster around formula 
(36). There is also good agreement at supersonic 
jet velocities. As noted before, the NASA Lang- 
ley data 10 is consistently 1 to 2 dB higher than the 
Tanna et a/. 33 data. The calculated results appear 
to lie right between these two sets of measurements. 

6. Concluding Remarks 

In this work a semi-empirical theory for the pre- 
diction of the fine scale turbulence noise component 
of high speed jets is developed. The fine scale tur- 
bulence noise is the dominant noise component in 
the sideline and upstream directions. By comparison 
with experimental measurements over a wide range 
of jet velocity and temperature ratios, it is found 
that the theory can provide very accurate noise pre- 
diction. 

The present theory is semi-empirical. Two crit- 
ical components of the theory are the model two- 
point space-time noise source correlation function 
and the k — e turbulence model. Both models con- 
tain a fair amount of empiricism. There is, therefore, 
room for improvements. Theoretically, it would be 
desirable to be able to derive the two-point space- 
time correlation function from first principle and/or 
to replace the k — e model by a less empirical tur- 
bulence theory. But such an endeavor is beyond the 
scope of the present work. 

Despite the inherent semi-empirical nature of 
the present theory, the accuracy of the prediction 
over such a broad range of jet velocities and tem- 
peratures makes it superior to other methods in the 
literature and in practice. The present method can 
provide accurate spectral prediction that is vital to 
applications in which the use of the perceived noise 
level (PNdB) is required. The prediction of noise 
spectrum has proven to be most difficult in the past. 
Many existing theories are capable of calculating 
only noise intensity. From a practical standpoint, 
we believe the present work should be of interest to 
the aeroacoustics community in need of a better jet 
noise prediction methodology. 
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Figure 1. Random motion of gas molecules in the moving frame 
fixed to the mean velocity V. 



Figure 2. Random motion of blobs of fine scale turbulent fluid in 
a frame of reference moving with the mean velocity V. 



Time detay, r, /isec 

Figure 3. Two-point space-time correlation of the axial velocity 
component in a jet. — — —measured data, Davies et al. (1963), 
model function. 



Figure 4. Schematic diagram showing the difference in ray paths 
between source points Xi and X 2 from a far field point x with 
spherical coordinates ( R , ©, <f>). 



Figure 5a. Cold jet noise spectra at © = 60 deg., = 1.0. 

(a) Mj = 2.0, (b) M } = 1.49, (c) Mj = 0.9, (d)°° M ; = 0.7, 

(e) Mj — 0.5, (f) Mj = 0.3. experiments, Ref [10],[32] 

equation (36). 
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